Compact vortexlike solutions in generalized Born-Infeld model 
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We study vortex-like solutions in a generalized Born-Infeld model. The model is driven by two 
distinct parameters, one which deals with the Born-Infeld term, and the other, which controls the 
presence of high-order-power term in the covariant derivative of the Higgs field. We numerically 
solve the equations of motion and depict the main vortex features, for several values of the two 
parameters of the model. The results indicate the presence of compact vortex, when the parameter 
responsible for the high-order-power in the derivative increases to sufficiently large values. 
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I. INTRODUCTION 

This work deals with the presence of vortices in mod- 
els described by Abelian gauge field coupled to complex 
scalar field. In the standard case, the scalar field self in- 
teracts via the presence of the Higgs potential, capable 
of inducing spontaneous symmetry breaking. Moreover, 
the gauge and scalar dynamics and interactions are con- 
trolled by the Maxwell term, and by the usual covariant 
derivative squared. The system is well known, and capa- 
ble of supporting vortex solutions [l[ , which are of general 
interest to physics 0, Q . 

Recently, several authors have studied models with 
modified dynamics, sometimes known as k-field models. 
The basic idea is to change the dynamics, allowing for the 
presence of terms containing higher-order power on the 
derivative of the field. Motivations to investigate these 
models appear in cosmology attempting to con- 

tribute to understand the present accelerated expansion 
phase of the universe, and to study issues of strong in- 
teraction and other subjects of current interest to high 
energy physics (l^-[30|. 

An interesting point of such generalized models is that 
they may contain defect solutions with compact support 
[IJ-js^l, known as compactons. If topological defects are 
seen as field-theoretical realizations of fundamental parti- 
cles, then compactons are truly localized solutions in the 
sense that they are absent of the typical infinite solitonic 
tails, reaching their vacuum values at finite distance, as 
opposed to conventional solitons. As such infinite tails 
are more than a modeling nuisance, compactons could 
be able to model more realistically nature. Compactons 
first appeared in [3l| and have gained many interest in 
high energy physics and in other recent investigations 
[Hi- (Hi J including compact bosonic stars [36j . 

In a recent work, some of us studied the presence of 
compact vortex in the Maxwell-Higgs model, modified to 
allow for the presence of a term with the fourth-order 
power in the covariant derivative of the complex scalar 
field [st'I . The main result suggests that the vortex solu- 
tions exhibit compact features when the term introduced 
to modify the dynamics is important. Motivated by this 
result, in the present work we go further on the sub- 



ject, and we investigate rotationally symmetric solutions 
obtained from a specific k-ficld theory. This theory is ob- 
tained from the Born-Infeld-Higgs Electrodynamics, and 
it is now endowed with an extra kinecticlike contribution, 
which is introduced to allow for the presence of a term 
depending on higher-order power of the covariant deriva- 
tive of the scalar field. Born-Infeld (BI) theory, originally 
introduced to deal with the problem of the divergence of 
the electron self energy in Classical Electrodynamics [s^ 
is singled out among the class of nonlinear electrodynam- 
ics by its special properties under wave propagation [s^ 
and it currently receives a special attention as it arises 
in the low-energy limit of the string/D-Branes physics 
poj - 42 1 . Moreover in the defects context this action has 
been found to support global string solutions [33| . 

The aim of this paper is twofold. First to study how 
the extra kinecticlike term contributes to the compacti- 
fication of the vortexlike structure present in the model 
with gauge fields, as such solutions have been already 
found in some k-field field models [s^l as well as in theo- 
ries with V-shaped potentials (see e.g. And second 
to explore how the BI gauge field aff'ects the existence of 
solutions which, as shown below, cease to be defined be- 
low some critical value of the BI parameter /3 = /3c, as it 
has been also suggested to happen in other models pi ), 
indicating that a high enough amount of nonlinearity in 
the gauge field generically leads to a breakdown in the 
existence of the solutions within this kind of nonlinear 
actions. 

We organize the paper as follows: in the next Sec. |lTl 
we introduce the generalized Born-Infeld-Higgs (gBIH) 
model and make the basic considerations which allows to 
get to the equations of motion to describe the vortex con- 
figurations. In Sec. mil we perform the numerical study 
and depict the main features of the vortices. Finally, in 
Sec. IIVI we include comments and conclusions. 



II. EQUATIONS OF MOTION 

Here we deal with the gBIH model in (2, 1) dimensions. 
The Lagrange density has the form 



C^Cbi + K{X)~V{\^\) 



(1) 



2 



The (non-zero) diagonal elements of the metric are 
(1,-1,-1), and the three terms in the above Lagrange 
density are given by 



(2) 



We fix ^4" = as a proper gauge choice and use it from 
now on. 

As usual, we consider the rotationally symmetric 
Ansatz for the fields 



(j) (r, 6) = vg (r) e' 



(12) 



and K{X) is a function of X to be given below, where X 
stands for 



X = \D„ 



(3) 



and -Dp0 represents the covariant derivative of the scalar 
field 



A(r, 



■(a(r) -n) , 



(13) 



where r and 9 are polar coordinates and n is an integer, 
the winding number (vorticity) of the solution. In this 
case, the rotationally symmetric equations of motion are 



and V^(|0|) stands for the potential 



(4) 



(5) 



In the above expressions, we use e and A as real and pos- 
itive coupling constants, and v counts for the symmetry 
breaking parameter, real and positive. 

For simplicity, we introduce the mass scale M and 
make the scale transformations: x^^ /M^ (j> ^■ 

Ml/20, ^ M^^^A'', P M^/^P, e ^ A/i/^e, 
V — > M^/^u and A — > Af^/^A. In this case, we get 
X M^X, K APK and C M^Cg, and the 
above model is now described by the Lagrange density 
Cg, which has the same form ([T]), but now all the coordi- 
nates, fields and parameters are dimensionless quantities. 

The procedure to search for vortex solutions requires 
that we investigate the equations of motion. In the 
present case, they can be written in the form 



a, 



[iz' 



= r 



(6) 



d I \ da 
dr \ r dr 




-1/2N 



(14) 



r dr \ dr 



v'^K — — I ( — 
dr dr \\ dr 

where X has now the form 
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2 



? 9 

9 a 



dg 



X 



IDc 



-V 




(16) 



We make the investigation specific, choosing the fol- 
lowing forms of K{X) and t^(|0|) 



K{X) = X - aX'^ 



(17) 



(18) 



KxD^D^cl, + KxxdpXDy = ~% . 
Here we use Kx = dK/dX, Kxx = d^K/dX^, and 



J" = ie UD'^(t> - Kx 



1 



2/32 ^^"^ • 



The static Gauss Law can be written as 



(7) 



-2e2.4" \6\^ 



K 



X 



where i? = ^ x A, E 
now the form 



_^^0 g012 



7^ 



(8) 
(9) 

(10) 

1 and TZ has 
(11) 



Note that we have introduced a new factor in the po- 
tential, depending on a, the parameter which controls 
the extra term in K[X), responsible for the generalized 
modification of the model, as it appears above. The rea- 
son for this extra term is explained as follows: when the 
scalar part in the lagrangian ([T]) is considered, with (fT7|) 
for the kinetic term and ([S]) for the standard potential, 
the equations of motion arc written as 



(19) 



where the new potential U{\(j)\) can be shown to be given 
by M 



TJiu\\ 1 y/l + 12«nH) - 1 

a 



(20) 



1 x/i+to-i ^ For simphcity it 



which implies that C/(0) 

is desirable to make the constant U (0) independent of a 



3 



and choosing for example f/(0) = 1/2 this leads to the 
modification 



l + -a]Vm), 



(21) 



which brings the potential (|18l) . As shown in Ref. [37[ this 
choice is a good way to investigate the compact profile in 
theories without the BI term. We take a nonnegative, to 
avoid problem with the positive definiteness of the energy 
of field configurations, as the energy density (23) below 
suggests; see also Ref. [23]. 

In this case, for the specific choice above, the rotation- 
ally symmetric equations of motion become, after taking 
e = f = A = 1, for simplicity: 




(22) 



/ 



0.5 



FIG. 1: We take a = \ and /3 = 2, and we plot a(r) with solid 
(black) line and g(r) with dashed (red) line, as a function of 
r. We also plot the standard solution, for a = and = 0, 
with a(r) as dotted (black) line and gir) as dash-dotted (red) 
line, for comparison. 



The energy density must vanish asymptotically, as a con- 
dition to make the energy finite. So, we get 



l + 2a 



2a— — 

dr dr 







dr ) 






g'aA 




J.2 I 




(23) 



These are the equations we have to deal with in the nu- 
merical investigation. We recall that the limit /3 — >■ oo 
leads us back to the model studied in Ref. [37[. Also, if we 
further take the limit a — )■ 0, we get back to the standard 
Maxwell-Higgs model. We also note that if we include the 
Chern-Simons term in the model ([T]), the presence of a 
generalizes the model studied before in ^43] , searching for 
electrically charged vortex solutions in the BI model, in 
the presence of the Chern-Simons term. However, we will 
not take this route here, that is, we will not include the 
Chern-Simons term in this investigation. 

To study the presence of vortex solutions to equations 
and one needs to know how the fields a{r) and 
g{r) behave at the origin, and in the limit r — > oo. Near 
the origin they should obey 



g (r — > 0) -> and a (r ^ 0) 



1 



(24) 



where we fix n — 1, for simplicity. The asymptotic 
boundary conditions are obtained from total energy, 
which is given by 



(25) 



J e{r) d^r, 
where the energy density has the form 



lim 10 (r) 



lim |D(^| = 



(27) 



(28) 



From these conditions, using the rotationally symmetric 
Ansatz, we get, respectively, the asymptotic boundary 
conditions on the profile functions 



g{r go) 1 and a (r — > go) — > 



III. NUMERICAL SOLUTIONS 



(29) 



Let us now focus our attention to the vortex solutions. 
The equations of motion to be considered are ([^^ and 
((23)) . and g and a have to obey the conditions (l24l) and 
(ESI). 

We also plot solutions for the energy density, explicitly 
written as 



£ (r) = /3 




and for the magnetic field 



r dr 



(31) 



The strategy to numerically solve the equations of mo- 
tion, given the fact that the boundary conditions do not 
lie at the same point, is to employ a shooting method. 
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FIG. 2: We take a = 1 and /3 = 2, and we plot e{r) (upper 
panel) and B{r) (lower panel), as a function of r. Dashed 
lines show the same quantities, but in the standard case, for 
a = and = 0. 



FIG. 4: We take /3 = 6, and we plot e{r) (upper panel) and 
B{r) (lower panel), as a function of r, for a = 1 and 10, with 
dashed (blue), and dash-dotted (red) lines, respectively. 




0.5 1 1.5 2 2.5 3 



FIG. 3: We take /3 = 6, and we plot a{r) and g{r) as a function 
of r, for a — 1 and 10, with solid (black) and dashed (red) 
lines, and with dotted (black) and dash-dotted (red) lines, 
respectively. 



suitably adapted to the system of equations (f22|) - ([23l) . 
Shooting from the center we numerically determine the 
constants g'{0) and a"(0), up to some predetermined pre- 
cision, and which are of course a and /3-dependent. This 
strategy is similar to the one employed, for example, in 
Ref.[3l|. 

The model to be investigated has two parameters, a 
and /3, which we have to vary to investigate how the 
vortexlike solutions depend on them. The parameter (3 
controls the BI term, and for large values of /3, the BI 



term should behave like the Maxwell term. Also, for a 
very small, the modification coming from the dynamics 
of the scalar field should be insignificant; see Eq. ^T7\ . 

In order to illustrate the procedure, let us first consider 
the case a = 1 and (3 = 2. The results for a(r), g{r), 
e{r), and B{r) are all depicted in Figs. 1 and 2. We have 
investigated the case with a = 1/13 = 0.1 as well, and 
it can hardly be distinguished from the standard case, 
corresponding to a = 1/(3 = 0, which is also plotted in 
the same Figs. 1 and 2, for comparison. We note from 
these figures that the values a = 1 and (3 = 2 show 
practically no effect of compactification on the vortexlike 
solutions. 

Let us now concentrate in the case of a large (3, which 
we choose to he (3 = 6. In this case, we implement the 
numerical investigation with a comparison with the pre- 
vious work [37| . Thus, we firstly fix /3 = 6, and we vary 
a, taking the values a = 1 and 10, recalling that the 
degree of importance of the extra term in X increases 
with increasing a; see Eq. (|17p. The numerical results 
for the fields a{r) and g{r), the energy density e{r) and 
the magnetic field B{r) are shown in Figs. 3 and 4, and 
there we see that they behave as they do in the case 
where the BI term is changed by the standard Maxwell 
term; compare with results of Ref. [s^l- In particular, 
the results depicted in Figs. 3 and 4 clearly show the ef- 
fects of compactification of the vortexlike structure, as a 
increases. 

Since the numerical study is now under control, we 
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FIG. 5: We take a = 10, and we plot B{r) as a function of r, 
for /3 = 6, and 4, with dashed (blue), and dash-dotted (red) 
lines, respectively. 
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FIG. 6: We take a = 0, and we plot B{r) as a function of r, 
for /3 = 5 and 1, with dashed (blue), and dash-dotted (red) 
lines, respectively. 

have investigated several distinct cases numerically. Be- 
low we comment on some specific values of the param- 
eters a and /3, which we use to illustrate the general 
situation. 

It is of current interest to understand if the nonlinear- 
ities introduced by the BI term contribute to stress the 
effects of compactification of the vortexlike structure. To 
study this issue, let us fix a at a large value, say a = 10, 
for which we already know that the vortex tends to have 
compact profile. Then, we numerically investigate the 
model for /? at the values /? = 6 and 4. These results 
are shown in Fig. 5, where we only plot the magnetic 
field B{r), which we think illustrate correctly the issue 
under investigation. Figure 5 shows that the effects of 
compactification increases as /3 decreases, so the nonlin- 
earities introduced from the BI term for /3 decreasing to- 
ward the value /3 = 1 also contribute to compactification. 
To check if this effect depends on a, let us now investi- 
gate the case with a = 0, for some values of /?. These 



results are plotted in Fig. 6, for /3 = 5 and 1, and they 
show that the effects of compactification are still there, 
although they are now much less significant. 

It is worth discussing, as a conclusion extracted from 
the numerical calculations performed, the fact that the 
solutions, for a given a, cease to exist below some critical 
value /3c(a), i-e. they only exist in the range /? < /?c(q;) < 
oo. For example, if we choose the value a = 1 this crit- 
ical value is placed around /3c — 1-25 and for a = 10 in 
the range 3.5 < /3c < 4. This feature gets manifest, as (3 
approaches to /3c, in a quick increasing of the slope uq of 
a{r) near the origin, in order to satisfy the correspond- 
ing boundary condition at infinity, in such a way that at 
the critical value we have lim^_>^^ao(r ~ 0) — cxd and 
the computations break down. This is better seen in the 
evolution of the magnetic field, as it is defined as the 
derivative of the a{r) field, which grows without limit as 
we approach the critical value. This effect seems to be a 
generic feature of Born-Infeld actions for vortexlike solu- 
tions, as can be seen, for example, in [3] or [i^, and not 
something artificially introduced by the numerical pro- 
cedure employed, as the checking with other numerical 
recipes (e.g. relaxation methods) may show. 



IV. FINAL COMMENTS 

In this work we studied the presence of vortexlike so- 
lutions in a BI model coupled to scalar matter field de- 
scribed by generalized dynamics. The model is controlled 
by two distinct parameters, a, which describes the modi- 
fication in the dynamics of the matter field, and /3, which 
drives the BI generalization. 

We investigated several regions in parameter space, 
and we could conclude that for a given /3, the increas- 
ing of a increases the effects of compactification of the 
vortexlike solution. We could also see that for a fixed, 
the increasing of 1//3 (up to some maximum value l//3c) 
increases the effects of compactification, although in this 
last case, the effects of compactification are less signifi- 
cant. However, the general conclusion is that the increas- 
ing of a and 1//3 enlarges the effects of compactification, 
but compactification only occurs at very large values of 
a. Other quantities such as the energy density and the 
magnetic field are significantly affected during the evolu- 
tion with the BI parameter, an effect also observed, for 
example, in Ref. [32|. 
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